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Abstract The theory of diffusion in many-dimensional Hamiltonian system is ap- 
plied to asteroidal dynamics. The general formulations developed by Chirikov is ap- 
plied to the Nesvorny-Morbidelli analytic model of three-body (three-orbit) mean- 
motion resonances (Jupiter-Saturn-asteroid system). In particular, we investigate 
the diffusion along and across the separatrices of the (5, —2, —2) resonance of the 
(490) Veritas asteroidal family and their relationship to diffusion in semi-major 
axis and eccentricity. The estimations of diffusion were obtained using the Mel- 
nikov integral, a Hadjidemetriou-type sympletic map and numerical integrations 
for times up to 10 8 years. 

Keywords Chaotic motion, Chirikov theory, asteroid belt, Nesvorny-Morbidelli 
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1 Introduction 



The application of chaotic dynamics concepts to asteroidal dynamics led to the un- 
derstanding of the main structural characteristics of asteroids distribution within 
the solar system. It was verified that chaotic region are generally devoid of larger 
asteroids while, in contrast, regular regions exhibit a great number of them (see for 
instance, Berry 1978, Wisdom 1982, Dermott and Murray 1983, Hadjidemetriou 
and Ichtiaroglou 1984, Ferraz-Mello et al. 1997, Tsiganis et al. 2002b, Knezevic 
2004, Varvoglis 2004) It was soon accepted that chaos was related inevitably to 
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instability, which may be local or global. Subsequent investigations searched for 
initial conditions leading to instabilities in relatively short time. In many appli- 
cations, the determination of Lyapunov exponent on a grid of initial conditions 
was used to get quantitative informations on stability. The inverse of the largest 
Lyapunov exponent, called Lyapunov time, should be in some way linked to the 
characteristic time for the onset of chaos (Morbidelli and Froeschle, 1996). 

However, some investigations have shown that many asteroids exhibit inter- 
mediary behavior between chaos and regularity. The first registered case was the 
asteroid (522) Helga (Milani and Nobili, 1992). This asteroid was in chaotic or- 
bit with a Lyapunov time inch shorter than the age of the solar system, but it 
exhibited a long period stability. No significant evolution was observed in the or- 
bital elements of (522) Helga for times up to one thousand times its Lyapunov 
time. Since then, other asteroids have been shown to have Lyapunov times much 
shorter than the stability times unraveled by simulations (e.g., Trojans, cf. Milani 
1993). Currently, this behavior is known in literature as stable chaos (Milani et al. 
1997, Tsiganis et al. 2002a, Tsiganis et al. 2002b). Indeed, there is strong evidence 
that local instability does not mean chaotic diffusion, in the sense that nothing 
can be said about how much global or local integrals (or orbital elements) could 
change in a chaotic domain, even when a linear stability analysis shows rather 
short Lyapunov times (see Giordano and Cincotta, 2004, Cincotta and Giordano 
2008). 

Nesvorny and Morbidelli (1998, 1999) demonstrated that one source of stable 
chaos is related with three-body (three-orbit) mean-motion resonances (Jupiter- 
Saturn-asteroid system) . They observed that asteroids in these resonances exhibit 
a slow diffusion in eccentricity and inclination, but no diffusion in the semi-major 
axis. According to the estimates of Nesvorny and Morbidelli (1998), about 1500 
among the first numbered asteroids are affected by three-body mean-motion res- 
onances. 

The three-body mean-motion resonances are very narrow since they appear at 
second order in planetary masses, their typical width being ~ 10~ 3 AU, but they 
are much more dense (in phase space) than standard two-body mean-motion reso- 
nances of similar size. Nesvorny and Morbidelli (1999) developed a detailed model 
for the three-body mean-motion resonance and presented analytical and numerical 
evidence that most of them exhibit a highly chaotic dynamics (at moderate-to- low- 
eccentricities) which may be explained in terms of an overlap of their associated 
multiplets. By multiplet, we refer to all resonances for which the time-derivative 
of the resonant angle, o"p,pj,p s , satisfies 

&P,Pj,ps = rnjkj + m s \ s + mk + pui + pjuij + p s ^ s ~ 0, (1) 

for given (mj,ms,m) £ Z 3 /{0}. In ([I]) the A's and cc's denote, as usual, the mean 
longitudes and perihelion longitudes, respectively; {p,pj,ps) £ ^ 3 are integers 
such that '^2 i (mi + Pi) = for i ranging over three bodies (Jupiter, Saturn and 
the asteroid). 

We will be dealing in this paper with the case (mj, mg, m) = (5, —2, —2). This 
three-body resonance seems to dominate the dynamics of, for instance, the aster- 
oids (3460) Ashkova, (2039) Payne-Gaposchkin and (490) Veritas (see Nesvorny 
and Morbidelli 1999). In the case of the first two of those asteroids (with relatively 
large eccentricity, ~ 0.15 — 0.20), their behavior looks regular over comparatively 
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long time-scales (typically ~ 1 — 10 x 10 3 years) while in case of (490) Veritas 
(with eccentricity, ~ 0.06) its dynamics looks rather chaotic over similar time- 
scales. The determination of the age of (490) Veritas family has been the concern 
of some authors who studied stable chaos (Milani and Farinella 1994, Knezevic 
1999, Knezevic et al. 2002, Knezevic 2003, Knezevic et al. 2004, Tsiganis et al. 
2007, Knezevic 2007, Novakovic et al. 2009). 

Herein, we investigate chaotic diffusion along (and also across) the above men- 
tioned three-body mean-motion resonance by means of a classical diffusion ap- 
proach. We use (partially) the formulations given by Chirikov (1979). That for- 
mulation are developed to study specifically Arnold diffusion or some kind of 
diffusion that geometrically resembles it initially called Fast- Arnold diffusion by 
Chirikov and Vecheslavov (1989, 1993), as well as the so called modulational diffu- 
sion (Chirikov et. al, 1985). However, although from the purely mathematical point 
of view several restrictions should be imposed, there are many unsolved aspects 
regarding general phase space diffusion (see for instance Lochak 1999, Cincotta 
2002, Cincotta and Giordano 2008). 

The structure of the Hamiltonian used by Chirikov in his first formulation is 
similar to the Hamiltonians obtained with the perturbations theories of Celes- 
tial Mechanics. In particular, the Hamiltonians of analytic models of the three- 
body mean-motion resonances are directly adaptable, with some restrictions, to 
Chirikov's formulations. 

Let us mention that some progress has been done in the study of Arnold diffu- 
sion, particularly when applied to simple dynamical systems, like maps, the latest 
ones are for instance, the works of Guzzo et. al (2009a), (2009b), Lega (2009). 
However the link between strictly Arnold diffusion and general diffusion in phase 
space is still an open matter. Indeed, Arnold diffusion requires a rather small per- 
turbation, when the measure of the regular component of phase space is close to 
one. Thus, as far as we know, almost all investigations regarding Arnold diffu- 
sion involves relatively simple dynamical systems like quasi-integrable maps. In 
more real systems, like the one investigated in this paper, the scenario is much 
more complex in the sense that the domain of the three body resonance is almost 
completely chaotic. 

Finally, this work is justified by the fact that an application of all those theories 
to real astronomical models is still needed. 

In Sect. 2, we summarize the general problem of computation of the diffusion 
rate along the resonance and we discuss the limitations and difficulties to follow 
Chirikov approach in case of this particular three-body mean-motion resonance. 
Section 3 is devoted to the resonant Hamiltonian (given in Nesvorny and Morbidelli 
(1999)) and its application to the (5, —2, —2) resonance. In Sect. 4, we construct the 
simplified (or two-resonance) and complete (or three-resonance) numerical models 
used in our investigations. Moreover, informations about the algorithms and initial 
conditions used for numerical integrations and the procedure for estimation of the 
diffusion are also considered in this section. In Sect. 5, we discuss the numerical 
results on diffusion in the (5, —2, —2) resonance. In this application of the Chirikov 
theory, we are concerned with the role of the perturbing resonances in the diffusion 
across and along the (5, —2, —2) resonance and their relationship to diffusion in 
semi-major axis and eccentricity. Finally, in Sect. 6, we investigate the behavior 
of the asymptotic diffusion decreasing the intensity of the perturbations in the 
(5, —2, —2) resonance. In this case, we are interested in the study of the diffusion 
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under the action of an arbitrarily weak perturbation considering scenarios close to 
that of the Arnold diffusion. 

2 Chirikov's Diffusion Theory 

In this section we give Chirikov's (1979) as well as Cincotta's (2002) description 
of diffusion theory in phase space in order to provide a self-consistent presentation 
of the subject. Since most of the results and discussions given here are included in 
at least these two reviews, we just address the basic theoretical aspects. 

Let us consider a Hamiltonian system having several periodic perturbations 
that can create resonances. The initial conditions are chosen such that the system 
is in the domain of a main resonance, called guiding resonance. The term of pertur- 
bation corresponding to the guiding resonance is separated from the others, which 
will be called of perturbing resonances. The Hamiltonian has the following form 

H = H (I) + eV G (I) cos (m G • 6) + eV (I, 0) , (2) 

with 

eV = e J2 (I) cos (m • 6») , (3) 

where Vq and mg are, respectively, the amplitude and resonant vector of the 
guiding resonance, V m and m are, respectively, the amplitude and resonant vectors 
of perturbing resonances. Here (1, 0) are the usual TV-dimensional action-angle 
coordinates for the unperturbed Hamiltonian Hq (N > 3), the vectors mg, m £ 
Z JV /{0} and Vq, V m are real functions. The small parameter perturbation, e, is a 
real number such that t < 1. The resonance condition is fixed by 

5(r)=m G . W (I r ) = 0. (4) 

The surface S (V) = in the action space, is called resonant surface. 

2.1 Dynamics of the guiding resonance in the actions space 

Let us first consider the simple case of one single resonance, that is, let us assume 
that all Vm = for m ^ m G , and we chose initial conditions close to the separatrix 
of the guiding resonance. In the uj-space, the resonance condition mc-u r = has a 
very simple structure, just a (N — l)-dimensional plane the normal of which is the 
resonant vector m.Q. In the I-space, m G • w r = leads to the (N — l)-dimensional 
resonant surface S (V) =0, whose local normal at the point I = V is 

n r = (^[mc-wfl)]) (5) 

In addition, we consider the (N — l)-dimensional surface Ho (I) = E (in I- 
space) and, if we suppose that u> (I r ) is an one-to-one application, we can also 
write Ho{oj) = Ho (I (<*>)) = E (in uj-space). 

The manifolds defined by the intersection of both resonant and energy surfaces 
has, in general, dimension N — 2. By definition, the frequency vector oj is normal to 
the energy surface in I-space, since it is the I-gradient of Hq. The latter condition, 



Chirikov- Arnold in the (5, —2, —2) resonance 



5 



together with the resonance condition Eqn. ([3]), shows that the resonant vector 
m<3 lies on a plane tangent to the energy surface at I = V. Furthermore, the 
equations of motion (only with Hq and the guiding resonant term) show that I 
is parallel to the constant vector mg. Thus the motion under a single resonant 
perturbation lies on the tangent plane to the energy surface at the point I = V in 
the direction of the resonant vector. 



2.2 Local change of basis 

Now, let us introduce a canonical transformation (1,0) — > (p,tp) by means of a 
generating function 

N / N \ 
F(p,0)=J2( + Pkm Oi, (6) 

i=l V k=l / 

where fi^. is a N x N matrix with nn = (mg),. The transform action equations 
are 

N N 

h = li + ^ PkVki, i>k = Pkfii- ( 7 ) 
fe=i e-i 

The phases V>fc,fc = 1,---,N are supposed to be non degenerate, i.e., ^j 3 - ^ 0. 
As Cincotta (2002) has shown, this transformation should better be thought as 
a local change of basis rather than as a local change of coordinates. The action 
vector whose components are (ij — /J) in the original basis {uj,j = 1, . . . , iV}, has 
components pj in the new basis = 1, • • • ,N\ constructed taking advantage 

of the particular geometry of resonances in action space. 

We choose, fi 1 = mi = mg and since the vector mg is orthogonal to the 
frequency vector ui r (due to the resonance condition), it seems natural to take 
fi 2 = u r I |oj r | . The remaining vectors of the basis are fi k = e^, k = 3, . . . , N, the 
vectors are orthonormal to each other and to /i 2 . Let us define one of the e^, 
say e s , orthogonal to the normal n r to the guiding resonance surface. In general, 
all the vectors e k will be orthogonal also to fi 1 , except e s . In general, e s will not 
be orthogonal to mg. Then, considering N = 3 and since p = PifA^i = 1, . . . , 3, we 
can say that p\ measures the deviations of the actual motion from the resonant 
point across the guiding resonance surface, p3 measures the deviation from the 
resonant value along the guiding resonance, while p 2 measures the variations in 
the unperturbed energy. 

For N > 3 degrees of freedom the subspace of intersection of the two surfaces 
leads to a manifold of N— 2 dimensions. Following Chirikov (1979), this subspace is 
called diffusion manifold. The N — 2 vectors locally span (at the resonant value) 
a tangent plane to the diffusion manifold called the diffusion plane. Then, in the 
new basis, the action vector may be written as: p = pimg +p2<*> r / |a/ r | +q, where 
q is confined to the diffusion plane q = Qk e k with q k = p^ for k = 3, . . . , N . 

We write now the Hamiltonian ([23) in terms of the new components of the 
action. Expanding up to second order in p k , using the orthogonal properties of 
the new basis, recalling that tpi is the resonant phase and neglecting the constant 
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terms, we obtain for (fc, £) ^ (1, 1) 

2 N N 

k=ie=i u 



with 



1 N N dw r 
M^=EE"H^, (9) 

i=l.j=l J 

M- G = M^ =mG ^ mG ^ (10) 

where we have written Vq, V (i/>) instead of Vq (p) , V (p, ip). These functions are 
evaluated at the point I = V or p = 0. 

In absence of perturbation (V = 0), the components Pk,k = 2, . . . , N are inte- 
grals of motion, which we set equal to zero so that I r is a point of the orbit. Then 
the Hamiltonian (JSj) reduces to 

Hfa,if>)rsHi(pt,1>i) + eV(il>), (11) 

where 

#1 = 2^+6% cos Vi (12) 

is the resonant Hamiltonian associated to the guiding resonance. It is a simple 
pendulum. Note that the stable equilibrium point of the pendulum is ipi = ir if 
MqVq > 0, or iii = if M G V G < 0. 

To transform the phase variables, we take into account that the dot product is 
invariant under a change of basis. Recalling that ip/, = Yli^ke^e^ then if v denotes 
the vector m in the new basis, we have: tp m = m-9 = v ■ xp, where = y\ ^i^tk- 
As we can see, while the m k are integers, the quantities v k are, in general, non- 
integer numbers, due to the scaling of the phase variables. 



2.3 Changes due to perturbation 

As mentioned above, for V = the p k are integrals of motion and since Hi is also an 
integral, we have the full set of N unperturbed integrals: Hi,p2,qk,k = 3, . . . , N. 
But if we switch on the perturbation, these quantities will change with time. 
This can be seen using the equations of motion for the Hamiltonian (JSj) , where 
tpj = dH/dpj,j = 1,...,N. Performing derivatives and integrating, considering 
that for V = 0, pi (£ 2) are constants and pi = M G ipi — X^fc2 Tfl^Pt' we obtain 

N , v 

i >k (t) = \^\t8 2k + j:( 1 ^- J ^)p e t+^i(t)+^ k0 , k>i (is) 

1=2 V ' 

where Sij is the Kronecker's delta and tpjo is a constant. To get tp m (t), we evaluate 
the dot product J^i v i^>i 



tp m = m ■ 6 = V ■ l/> = £mtpl (t) + Wmt + Pm + K m , 



(14) 



Chirikov- Arnold in the (5, —2, —2) resonance 



7 



where 

JV 



Cm = V" v k (m) ^p-, uj m = m u> r , (15) 

AC— 1 

and /3 m is a constant and 



t=2 

The second relation of (|12p is obtained taking into account that 



m ' wT = ^X^i( m )Mij ' (/*2 | w 1) 



(17) 



and the fact that, since fi 2 is orthogonal to all fi i ,i ^ 2 and p2 ■ M2 = 1, the dot 
product only contributes to z = 2. Then, 

tj m = m • ijj 1 = U2 (m) |aj r J . (18) 

We are now ready to compute the time variation of the unperturbed integrals. 
From (|11[) and ([3|), for pf. = —dH/dipk^k ^ 1, we easily find 

Pk (t) « e ^ (m) ^ m sin(y3m (t) . (19) 

where V m = V m (V) . This equation holds for every component of the momentum 
p, except for p\. Since pi is not an integral, we use Hi, instead of pi. 

Chirikov (1979) calculated the total variation of Hi with the aim of construct- 
ing a whisker map to describe the Arnold diffusion. However, instead of it, we 
prefer, in the study of three-body resonances, to compute the evolution of the 
components of the momentum p by means of numerical integrations or, alterna- 
tively, by mean of a Hadjidemetriou-type sympletic map (see Sect. 4). However, 
we use a variation of Chirikov's construction to obtain a theoretical estimate of the 
slow diffusion. We then proceed and compute the total variation of p^. For details 
about the construction of the whisker map we refer to Chirikov (1979, Sect. 7.3) 
(see also Cincotta 2002 and the Appendix B of Ferraz-Mello 2007). 

If e is small enough, the phase space domains associated with all resonances 
present in ([3]) do not overlap. Then a standard procedure is to replace ipi (t) 
and ipi (t) by the values on the unperturbed separatrix and to solve analytically 
(|19[) . We make first the integration of (|19[) over a complete trajectory inside the 
stochastic layer assuming that ipi = tp{ x and K m = 0. Indeed, as mentioned 
previously for V = the pg(£ ^ 1) are integrals of motion and the phases tp m 
can be estimated considering pg (I ^ 1) = 0, such that K m = 0. Then, the total 
variations of the pj/s are given by 

Apk (t) w e Yl v k ( m ) v ™ I sin^m (t) dt, (20) 

where j)„ (t) = £m4>i x (t) + o; m t + /3 m . The estimate of integral into (f20|) is done 
considering the known solutions for the phase ipf x (t) obtained near both branches 
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of unperturbed separatrix of the pendulum H\. More details about these calcu- 
lations are given in the appendix of this paper. Here we only described the main 
steps and the final result for Ap k (t). 

Chirikov shown that the contributions of integral in (|20[) in both branches of 
separatrix are described in terms of the Melnikov integral with arguments 

(21) 

where the double sign indicates the both separatrix branches and Qq is the proper 
frequency of the pendulum Hamiltonian H\. In order to simplify the calculations, 
Chirikov considered only even perturbing resonances and the contribution of Mel- 
nikov integral with negative argument was neglected under the condition | A m | 3> 1. 
In contrast, the perturbations in the three-body mean-motion resonance model are 
non even and the arguments are small. Moreover, the asymmetry in the Nesvorny- 
Morbidelli model implies that the time of permanence of the motion near each 
separatrix is different. Thus, we introduce the factor Rj- which takes into account 
the difference in the time of permanence of the motion in each separatrix branch. 
Hence, after some algebraic manipulations the Eqn. (|20|) is rewritten as 

Apk ~7te i>fc (m) Q m sin (22) 

with 

Qm = VL [RtA 2{U{ (|Am|) + (1 - R T ) 4 2 | Cm | (- |Am|)] , (23) 

where ip^ = ip S m (t = t°) with ipf* {t = t°) = n. Equation (f22|) is a theoretical 
estimate for the total variation of the momenta p^'s inside the stochastic layer 
around of separatrix of the pendulum Hamiltonian Hi, and it is valid for non-even 
perturbation and for small A m . Estimations of the Melnikov integral, A 2 \^ m \ (|A m |), 
in terms of ordinary function can be obtained from the values of |A m | and |£m|- 
On the other hand, the factor Rx can be estimated from numerical experiments. 



±|A r 



2.4 The diffusion rate 

In Chirikov's theory of slow diffusion, each resonance has a role in the dynamics 
of system. The main resonance, that is the guiding resonance, defines the domain 
where diffusion occurs. The stronger perturbing resonance is called layer reso- 
nance. That resonance perturbs the guiding resonance separatrix and it generates 
the stochastic layer and its properties (width, KS-entropy, etc.). Thus, the layer 
resonance controls the dynamics across the stochastic layer. The weaker perturb- 
ing resonances are called driving resonances. They perturb the stochastic layer and 
control the dynamics along the stochastic layer. Then, the driving resonances are 
responsible for the drift along the stochastic layer, i.e., the slow diffusion. We are 
interested in obtaining an analytical estimate for the slow diffusion. To fulfill this 
task, we will estimate the diffusion in the actions whose direction is given along 
the stochastic layer. 

We introduced the slow diffusion tensor 

% = ^P i,J = S,...N, (24) 
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where T a = ln(32e/w s ) / Oq is the characteristic time of the motion within the 
stochastic layer of the guiding resonance (equal to half the period of libration or 
to one period of circulation of ipi near the separatrix) and the average in the nu- 
merator is done over successive values of <p m . Here w s is the width of the stochastic 
layer given by 

Ws = -T^Wmi — - — £ Qm L > 0. (25) 

(see Sects. 6.2 and 7.3 of Chirikov 1979). In the last equation, the subscript L 
indicate the layer resonance. The components of the diffusion tensor (|24p are es- 
timated using the Eqn. ([22)1 . Hence, because of dependence with the phase <p mD , 
the average in (|24[) depends: (1) of the correlation between successive values (p mD 
when the system approaches the edges of the layer; (2) of the possible interferences 
of several driving resonances. However, the analysis done by Chirikov shown that 
the terms that contribute to the diffusion must have the same phase <p mD (see 
Sect. 7.5 of Chirikov 1979 and Cincotta 2002 for more details). Hence, using (|22[) 
the diffusion tensor components in (|24[) are described as 



Dij = ^ 2 J~] n (mo) uj (m D ) Q mD s'm 2 <p mD . (26) 

Terms with different mr> are averaged out. 

Now, there still remains the problem of estimating sin 2 <p mr> . To solve this 
problem we need to consider that the structure of the stochastic layer affects the 
motion of the system. In fact, studies of the slow diffusion theories have shown that 
the stochastic layer is formed by two different regions. The first, more central, near 
the unperturbed separatrix, is totally chaotic. The second, more external, near the 
edge of the stochastic layer, includes domains of regular motion forming stability 
islands. When the solution approaches the edge of the stochastic layer, it could 
remain rather close to the neighborhood of those stability islands for long times. 
This phenomenon, called stickiness, leads to a reduction in the diffusion rate (for 
more details about the stickiness phenomenon see the recent work of Sun and Zhou 
2009 and references therein). Thus, near stability islands some correlations in the 
phases arise, which dominate the motion across and along the stochastic layer. In 
this case, the evolution of phases (p mL and ip mD cannot be random simultaneously, 
and their correlation decreases the diffusion rate (see Chirikov 1979, Cincotta 
2002). 

In order to estimate the correlation between sin 2 <p mr> and sin 2 <p mL , we use the 
so called reduced stochasticity approximation, introduced by Chirikov (1979) like an 
additional hypothesis. Hence, the theoretical rate of diffusion given by (|26[) may 
be now evaluated and has the form 

D n = on^rr /J Rm ° v% ( m£) ) v i Qm D i, j = 3, . . . , N. (27) 

The Eqn. (|27p is an estimate for the theoretical diffusion inside the stochastic 
layer. The diffusion coefficient includes two parameters that reduce the diffusion 
rate: itydue to non-even perturbations and R mD due to the reduced stochasticity 
approximation. The expression given here for the diffusion tensor is different of 
that given by Chirikov because of the introduction of the parameter Rj- and by the 
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possibility of having a small argument in the Melnikov integral. Moreover, we have 
considered that the reduction factor due the reduced stochasticity approximation 
is different for each driving resonance, while Chirikov considers the same value for 
all of them. 



3 Application to 3-body mean-motion resonance 

The Hamiltonian, in the extended phase space, associated to a given (mj,mg,m) 
resonance, in Delaunay action-angles variables, is 

H = ~ 21? + njAj + nsAs + VjITj + VsIIs + Vsec + Vres ' (28) 

where 

X, m, Xj, wj, Xs, ra?s (29) 

are the mean longitudes and longitudes of the perihelions of the asteroid, Jupiter 
and Saturn, respectively, and 

L = y/E-n = y/a{yi - e 2 - l) ,Aj,IIj,As,n s (30) 

are the actions conjugated to them. The frequencies nj,vj,ng,vg are the mean- 
motion and perihelion motions of Jupiter and Saturn, respectively. 

The first term in ([28} describes the Keplerian motion of the asteroid and the 
terms proportional to the planetary actions extend the phase space to incorporate 
the motion of the angles A, zu, Xj, zuj, Xg, tag in the unperturbed Hamiltonian. 
Details concerning the derivation of this Hamiltonian are given by Nesvorny and 
Morbidelli (1999), whose main results and formula are used in this paper. Note 
that this Hamiltonian does not satisfy the convexity condition, however, this fact 
should not be a restriction for the application of Chirikov's diffusion theory. 

The perturbing function, following Nesvorny and Morbidelli (1999), has been 
splitted into its secular and resonant parts 

Vsec = — Psec ("res) <^ Cj J eg S COS (ijVJj + is^S + (31) 

kj,k s ,k,ij,is,i 

Pres = — Yl Pres ( Qres ) eke j' e S S COS i^PJ'PS ) (32) 

kj,ks,k,pj,ps,P 

where, a res = fflres/aj is the semi-major axis corresponding to the exact resonance, 
a p,pj,ps = m jXj + mgXs + mX+pvj+pjvjj +pszug, fij is Jupiter's mass, e, ej, e$ 
are the asteroid, Jupiter and Saturn's eccentricities, respectively, and P sec («res), 
-Pres (ares) are given functions that are linear in Saturn's mass (see bellow). The 
harmonic coefficients satisfy d'Alembert rules, ij+ig+i = 0, mj+ms-\-m+p+pj+ 
PS = and the series are truncated at some order in \kj\ + \kg\ + \k\, + + \i\ 
and \mj\ + \mg\ + \m\ + \p\ + \pj\ + \ps\- Next, we reduce the secular part (|3ip to the 
quadratic term in asteroid's eccentricity in order to break the degeneracy of the 
unperturbed Hamiltonian, and introduce in (|28[) the new action-angle variables, 
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defined by 



and 



I' = (JV, Nj, N S , n, II j, n s ) (actions) (33) 

0' = {v,vj,v s ,w,wj,w s ) (angles) (34) 

= m j\ j + m s \s + m\, uj = \j, v s = A s , (35) 

mN, Aj = mjN + Nj, A s =m s N + N s . (36) 



The variables {IIj,vjj,ng,vjs) remain unchanged. We recall that the resonant 
perturbation (|32[) does not depend on vj and v$ (so that Nj,N$ are constant 
that we can take as equal to zero). Let us write 

\ = {N,n,n,,,n s ), e = {y,w,wj,w s ). (37) 

Eliminating the constant terms, the Hamiltonian (|'28p may be written 

H(I,0)=H o (I) + V(I,0), (38) 

where 

i „ f. n ^ 2 



H ° ( : ) = - 2m2A r2 " A) ( 1 + ^ J + (mjnj + m s n s ) N + vjllj + u s n s , (39) 
is the unperturbed Hamiltonian and the perturbation is described by 

V(I,0) = ^/? m (I)cos(m-6»), (40) 

m 

with 

m= (1,P,Pj,Ps), (41) 

and 

/3m (I) = — J2 p res{ar es ) e k e k /e k s s . (42) 
k j , 5 , k 

Nesvorny developed a procedure allowing to obtain the coefficients (|42[) in terms 
of power series of the asteroid eccentricity only. In the last column of Table [1] are 
the coefficients calculated by Nesvorny for the guiding (G), layer (L) and driving 
(D) resonances used in our numerical experiments. 

We have considered the guiding resonance, defined by the vector mg = (1, — 1, 0, 0), 
the layer resonance, defined by the vector = (1,0,-1,0) and the driving reso- 
nance defined by vector = (1,0,0,-1). The unperturbed separatrices of those 
resonances in the plane a — e are shown in Fig. ([1) . 

The next step is to introduce the Chirikov variables (p,ip) allowing to have 
a separate representation of the actions across and along the resonance within 
the stochastic domain of the guiding resonance. The canonical transformation, is 
performed by the generating function ([6]), with a transformation matrix, fi, given 

by 
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Table 1 Old and new resonant vectors, and coefficients of the guiding (G), layer (L) and 
driving (D) resonances (The coefficients were taken from Nesvorny and Morbidelli, 1999). 





vectors m 


vectors v 


coefficients /3 m (xlO 8 ) 


G 


(1,-1,0,0) 


(1,0,0,0) 


45.59e - 32.24e 3 


L 


(1,0,-1,0) 


(0.55,0.66,0.76,0.70) 


-2.76 + 0.93e 2 


D 


(1,0,0, -1) 


(0.68,0.47,0.92,-0.70) 


1.18 - 0.38e 2 



1 


-1 








T" 


7" 


vj 






\u> r \ 


|w r | 


lu,rl 


2v s n r 2 u}l 




vjv S (n r 2 -nl) 


|v| 2 K- 


Isfl 


|q' r l 


kfl 


Iql 


V2vj 


V2vj 


V2W2 





2 Ivl 


2 Ivl 


V 



(43) 



where 

|q r | = + M 4 ) K - + Inl 2 , 

|n r | = ^K) 2 + K) 2 , |v| = ^J + 2K) 2 . 

Once the matrix of the transformation is defined, the new variables (p,i/') can 
be rapidly obtained using the relations (J7J) . In the new basis the arguments of the 
periodic terms change. The new vectors v defined by m • = v ■ ip are shown in 
Table [TJ in addition to the resonant vectors m and their respective coefficients. 

The procedure of previous section was applied in the Nesvorny-Morbidelli 
model, and leads to the Hamiltonian (JSJ) with iV = 4. The three perturbation 
coefficients /3q, 0l and /3rj of the guiding, layer and driving resonances, respec- 
tively, are calculated at the resonant values I r , which satisfies In the plane 
Nil the resonant condition ^ leads to a curve satisfying to 



n r + dn r + C 2 = 0, 



(44) 



where C\ {N r ) and C 2 (N r ) are given in terms of N r . Then, the solutions of 
(|44[1 can be obtained analytically for a fixed value of N r . However, in Nesvorny- 
Morbidelli model the coefficients /3 m 's are given as functions of the asteroid ec- 
centricity (see Table [T|). Therefore, we must use the definitions of Delaunay vari- 
ables (|30[) to determinate (a r , e r ). The resonant eccentricity is determined through 

e r = Jl - (I + n r /N r ) 2 , where (N r ,I7 r ) satisfies ©. The resonant semi-major 

axis is determinate using a r = (N r /2) 2 . 



4 Numerical Experiments 

In this section, we describe the numerical experiments done to investigate the 
diffusion across and along the stochastic layer of the three-body mean-motion 
resonance (mj, 7715, m) = (5,-2,-2) and its relations with the diffusion in semi- 
major axis and eccentricity. In these investigations the diffusion across will be 
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e o.i5 



0.10 



0.05 



3.170 3.172 3.174 3.176 3.178 3.180 

a 

Fig. 1 Unperturbed separatrices of guiding, layer and driving resonances in the plane (a, e). 
(Nesvorny and Morbidelli, 1999) 

described by the actions (pi,P2) and the diffusion along by the actions (ps,p4). In 
order to determine the time evolution of each action p, we use the equations of 
motion obtained from Hamiltonian ([8]) with N = 4. Then, for each value p/-(f), we 
use the equations of transformation ([7|) to obtain the respective values of N(t) and 
n(t) and the definition of the Delaunay actions in (|30p to obtain a(t) and e(t). 

Two main models were considered in the numerical experiments: (i) simplified 
(or two-resonance model) and (ii) complete (or three-resonance model) . In the first 
one, only one term of the perturbation - the layer resonance - is considered. In the 
complete model, two terms are considered: the layer and one driving resonance. 
In both cases the guiding resonance is given by mg = (1,-1,0,0). 

Two different techniques were used to construct the solutions. In a first set 
of experiments, the equations of motion of the Hamiltonian (|8]) were numerically 
integrated using the Burlish-Stoer method, for times in the interval 10 2 < ti n t < 
10 years. The results of these simulations were sampled with an output time 
step of 10 years. The simulations were done for the eccentricities 0.05 and 0.25 
with the initial conditions given on the separatrix of the guiding resonance (pro = 
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2^J\MqP^,P2q = 0,p3o = 0,p4o = 0,tp = 0). The main goal in these experiments 
was the study of the variation of the rate of diffusion across and along as a function 
of the total time of the simulations. Moreover, we investigate the correlations 
between the diffusion in the Chirikov actions p and the diffusion in semi-axis 
major and eccentricity. 

In the other set of experiments, the simulations were done using an Hadjidemetriou- 
type sympletic mapping (Hadjidemetriou 1986, 1988, 1991, 1993; Ferraz-Mello 
1997; Roig and Ferraz-Mello 1999, Lhotka 2009) defined by the canonical trans- 
formation (p n ,i/> TI ) — > (p n+1 , ) whose generating function is given by 

S (p n+1 ,i, n ) =J2i>?P? +1 +VH (p n+1 ^ n ) , (45) 



where r\ is the mapping step and the Hamiltonian is given by ([8]). The mapping 
equations are 

n+1 n dH (p" +1 ,V>") 
Pi = Pi~ V -g^n " ( 46 ) 



ipi =4>i +V i = l,2,3. (47) 



The procedure to determinate the semi-major axis and eccentricity for each point 
(p™ +1 , °f the trajectory is analogous to that discussed above. The goal of 

these experiments is to obtain the diffusion contour plots in the region of the 
5,-2,-2 resonance in the plane (a, e) (that plane is shown in Fig. [I]) for the two 
models considered. In this case, the total time of integration used is the 10 s years 
with i] = 10 years. The initial conditions are defined by the knots of a grid in the 
plane (a, e), on the rectangle (3.17 < a < 3.18)U.A., (0.01 < e < 0.30). The initial 
condition of the state vector po, for each point of the grid, was obtained using the 
transformation equations Q an d the definitions of Delaunay variables. The initial 
condition for the phases is tp^Q = 0,k = 1,...4. The use of the Hadjidemetriou 
map was instrumental allowing the computation of the solutions starting on each 
point of the grid which, otherwise, would demand an excessively large amount of 
CPU-time. The comparison of results provided by the map with those obtained 
by integrating the Hamiltonian flow, do not show significant differences in the 
numerical computation of the diffusion coefficient (see below) , at least for the two 
values of the eccentricity used (0.05 and 0.25). 

Finally, we need a numerical procedure to estimate the diffusion coefficient of 
each element of the set (pi,p2,P3,P4, a, e). In his investigations, Chirikov (1979, et 
al. 1979, 1985) used a particular method to determine the diffusion coefficient of the 
total energy H of the system. After Chirikov (1979), this procedure allows the pro- 
cesses that are really stochastic to be separated from those associated to bounded 
oscillations of periodic nature. Chirikov's procedure for experimental determina- 
tion of the diffusion coefficient consist in dividing the total time of simulation ti n t 
in N k sub-intervals of length (At) k and the calculation of the mean value, pi, for 
every sub-interval. The contribution to the diffusion rate for a given pair p im , sep- 
arated by interval of time (m — £) (At) k , is given by {pi m — Pi e ) 2 / \m — i\ {At) k . To 
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obtain the rate of diffusion, the contributions of the considered pairs are averaged 
over all the combinations m ^ I. That is, 

D k = 2 V fe m -P*t) 2 (as,-) 

The sub-intervals, used to estimate the mean values of quantities p.; , were obtained 
with k = 10 and the length (<4i)io = t int /10. 

The same procedure was used to determine the diffusion of the semi-major axis, 
Da, and of the eccentricity, D\. We have also estimated the eccentricity variation 
in these experiments using a definition of diffusion rate of the random walking 
type (see for example Eqn. (1241) 1: 

8e ~ \fbji~ t . (49) 



5 Results and discussion 

In this sections, we discuss the results obtained in the numerical experiments 
described above. In the discussion we will call action across to (pi,P2), and across 
diffusion to ^D\,Di), where we suppressed the superscript k. In the same way we 
call action along to (p3,Pi) and along diffusion to (1)3,1)4). 



5.1 The role of number of the perturbing resonances in the diffusion 

In his theory, Chirikov showed that the number of perturbing resonances is im- 
portant for the dynamics of systems with many-dimensional Hamiltonians. The 
results, in this case, repeat what is know from the general theory of Hamiltonian 
systems. In a system with two degrees of freedom, the resonances may be iso- 
lated by KAM tori, but for (N > 3) the dimensionality may allows, in principle, a 
solution to visit the whole phase space when t — > 00. 

Several experiments, using the Burlish-Stoer integrator, were done see the way 
in which the number of perturbing resonances in the diffusion behavior. Figure [2] 
shows the results for the diffusion coefficients Di, i = 1, 2, 3, 4 in the simplified and 
complete models as function of total integration time for eccentricities equal to 0.05 
and 0.25, . In the plots of Fig. [2] we see that the estimated diffusion increases in 
the low eccentricities up to a maximum reached for 10 3 — 10 4 years. This behavior 
is explained by the fact that the solution needs to fill the stochastic domain in the 
direction across to it. After that maximum, in the simplified model the diffusion 
coefficients for all actions decrease continuously. This decrease indicates that the 
variation of the momenta in both directions, across and along the stochastic layer 
are bounded (as the total time increase, only the denominator of (|48p grows making 
the result to decrease). As predicted by Chirikov's theory of slow diffusion, the 
actions p$ and p4 along the resonance do not evolve, notwithstanding the absence of 
topological barriers for its evolution. Without a driving resonance, there is no long- 
period evolution of the solution along the stochastic domain. In our experiments, 
a very distinctive reduction in the diffusion is observed in the case e = 0.25 after 
tint ~ 10 7 years. This behavior is likely due to a sticking of the solution to some 
regular domain. 
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The behavior of the diffusion in the complete model is more complicated. In 
the experiment with e = 0.05, the diffusion coefficients for the actions across the 
stochastic domain after 10 8 years are smaller than for the actions along it. This 
difference reaches approximately four orders of magnitude in this case and is due, 
probably, to the limitation of the motion across the stochastic layer imposed by 
its width. For e = 0.25 (right plot of Fig. [2]), the diffusion coefficients in the two 
models present almost the same characteristics observed for e = 0.05, except by 
the fact that, now, the diffusion in the actions along the stochastic layer, present 
a slow reduction with the integration time after 10 4 years. This behavior is likely 
due to the absence of overlapping of resonance at high eccentricities, in contrast 
with the case of low eccentricities, where the three resonances overlap (see Fig.[TJ). 



5.2 Diffusion in semi-major axis and eccentricity in the complete model 

The study of the previous section was completed with the computation of the 
diffusion coefficients for the orbital elements: semi-major axis and eccentricity in 
the complete model. Figure[3]presents the results. The results for the actions shown 
in this figure are the same shown in Fig. [2 but with a magnified scale. We see that, 
for large total times, there exist a correspondence between the diffusion coefficients 
of the actions across (jpi,pi) and of the semi-major axis, and between the diffusion 
coefficients of the actions along the resonance (j>3,pi) and of the eccentricity. This 



Chirikov- Arnold in the (5, —2, —2) resonance 



17 



behavior can be understood observing the geometry of resonance (5, —2, —2) shown 
in the Fig.[T] The separatrices of resonances are straight lines and the motion, along 
one of these separatrices, has constant semi-major axis and variable eccentricity. 
Following the discussion presented in Sect. 15.11 and the comparison done in the 
previous section for the simplified and completed models, we know that the drift 
along the separatrices only occurs if there is at least one driving resonance. Hence, 
the eccentricity diffusion is due to the driving resonance. 

As a complement to the previous discussion, we note that the variations in 
semi-major axis occurs in the horizontal direction, the same direction of the actions 
(pi,P2). The behavior of the diffusion in semi-major axis is similar to the diffusion 
of the actions across the resonance (pi,P2) and is bounded by the width of the 
stochastic domain. A consequence of this fact is that the diffusion coefficient in the 
semi-major axis is smaller than that for the eccentricity (in the complete model, 
the diffusion along is not bounded). 

The Fig. [3] shows the variation of the eccentricity calculated using initial con- 
ditions forming a grid in the plane (a, e) (the same grid of Fig. [5]) in t int = 10 s 
years. The results for the simplified model are shown in Fig. HJa). In this case, 
the larger variations in eccentricity occur for small values of the eccentricity. Two 
shallow maximums are formed, which are likely related with the eccentricity value 
at the intersection of the separatrices of the guiding and layer resonances (the only 
secondary resonance considered in the simplified model) . The results for the com- 
plete model are shown in Fig.[3Jb). In this case, the eccentricity variation reaches 
high values in the domain of low eccentricities - between O.Oland 0.125 - with a 
maximum for < e >~ 0.05. This maximum is certainly a result of the overlapping 
of the resonances in low eccentricities, forcing the actions along the resonance. 

In this model for mean eccentricities between 0.125 and 0.20 the variations 
are of the same order. The distributions observed in the Fig.[DJb) is in agreement 
with Nesvorny's unpublished data for 45 numbered asteroids of the (5,-2,-2) 
resonance (see the Table 2 in Nesvorny and Morbidelli 1998). The use of the 
models with only one perturbing resonance does not allow to get the distribution 
of the eccentricity variation observed in Fig. Hfb). 

5.3 The stochastic domain in the plane (a, e). Dependence on the initial 
conditions 

The diffusion coefficients were calculated on a large set of initial conditions to 
assess the domain where the solutions present stochastic behavior. The analysis 
was done using simulations over t int ~ 10 8 years, on the points of a grid of initial 
conditions in plane (a, e). A Hadjidemetriou-type sympletic mapping was used 
instead of expensive numerical integration to allow a large number of simulations. 
Figure [5] shows the contour plots of the diffusion coefficients of P4. It shows the 
stochastic domain of the guiding resonance (the light gray areas in Fig. [5]) . Note 
that the stochastic domain follows the geometry of the unperturbed separatrix of 
Fig. HJ Also note that the results for the complete model show a stochastic domain 
(a, e) larger than that observed for the simplified model. 

These differences are easily understood if we note the overlapping of the three 
resonances in the considered range of eccentricities. Figure [T] shows that the sepa- 
ratrices of the layer and driving resonances are, for almost all eccentricities, interior 
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e = 0.05 e = 0.25 




Fig. 3 Diffusion Coefficient for actions across and along, semi-major axis and eccentricity in 
experiments, obtained for complete model, for times from 10 2 up to 10 s years. Each point is one 
experiment with initial conditions upon the unperturbed separatrix of the guiding resonance. 



to the domain of the guiding resonance. At low eccentricities, however, the sep- 
aratrices cross one another. Thus, in low eccentricities, one solution crossing the 
chaotic neighborhood of the separatrix of the guiding resonance, also cross the sep- 
aratrices of the layer and driving resonances. The driving resonance acts pushing 
the actions along the guiding resonance. The magnitude of the push is determined 
by the phase ipmr> an d amplitude /?m D ■ 

At variance with the complete model, the simplified model presents very low 
diffusion, in low eccentricities, as seen in Fig. [2] In this case, the absence of the 
driving resonance (only the guiding and layer are considered in the simplified 
model) implies in the absence of evolution along the guiding resonance. 

A remarkable feature in both results is the formation of a wide region, in 
the central part of the domain of the guiding resonance, where the diffusion is 
negligible. The motion appears regular for initial conditions inside that region 
even when considering very long time spans. This result confirms what Nesvorny 
and Morbidelli (f 999) observed in surface of sections for eccentricity 0.20 using this 
same analytic model reduced to two degrees of freedom and two resonances. This 
is different from the situation observed in low eccentricities, where the separatrices 
of the resonances overlap. 
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Fig. 4 Variation of the eccentricity versus mean eccentricity for (a) simplified model and (b) 
complete model on a net of points in the plane (a, e) . 



6 Asymptotic behavior 



Chirikov theory of slow diffusion was constructed to study the diffusion under the 
action of an arbitrarily weak perturbations, and the diffusion coefficient was com- 
puted there using the asymptotic estimate of Melnikov's Integral. The asymptotic 
behavior of the three-body resonance model of Nesvorny and Morbidelli was stud- 
ied using the same technique devised by Chirikov. Figure [6] shows the variation 
of the diffusion coefficients along the resonance, for two different initial eccentric- 
ities (0.05 and 0.2), as functions of the parameter A mD = uJ mD /^G appearing as 
argument of the Melnikov integrals in Sect. 2.3 in the case the driving resonance 
m = (1, 0, 0,-1). Small values of Xm D ar e obtained decreasing the intensity of the 
guiding and perturbing resonances. 

The Hadjidemetriou-like mapping was used to allow us to compute the solu- 
tions over 10 10 years for a great deal of different conditions. The diffusion coefficient 
was calculated for initial conditions over the separatrix of the guiding resonance. 
A background value D^, to be used as reference, was also obtained with initial 
conditions in the central part of the guiding resonance, far from the separatrices. 
For small values of the perturbation, the motion in the central part of the res- 
onance domain is regular and the background diffusion appear as smaller than 
the diffusion shown by solution starting on the separatrices. For high values of 
the perturbation intensity, the motion is chaotic over the whole domain and the 
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Simplified Model 




3,1700 3,1720 3,1740 3,1760 3,1780 3,1800 



Fig. 5 Diffusion coefficients for action along j>4 for initial conditions in the interval 3.17 < 
ao < 3.18 U.A. for semi-major axis and 0.01 < eo < 0.30 for eccentricity for both simplified 
and complete models. The results were obtained for total integration time equal to 10 s years. 
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diffusion coefficients in the central part are not different form those of solutions 
starting on the separatrices. 

Figure [6] shows the diffusion coefficient D3 for e = 0.2 and for the low eccentric- 
ity case e = 0.05. The figures for the coefficient D4 are not shown since they are 
almost identical to those shown for D3. Figure [6] shows the 3 different possibilities. 

1. The first section of the figures, corresponding roughly to A mD < 2, is character- 
ized by complete chaos. For the smallest Xm D > one sees the same phenomenon 
discussed in Sect. [5~T1 the background diffusion appear small for some shorter 
runs because they do not cover the time necessary to allow the solution to fill 
the chaotic layer; but when time span grows, the diffusion values increase as 
expected. In this section, in general the diffusion coefficients for solutions start- 
ing in the central part or on the separatrices are equal showing that the whole 
resonance domain is chaotic. A few exceptions appear, as shown in Fig. Eft). 
In addition we may see in this figure, for A mD ~ 1, a sudden decrease of the 
background diffusion indicating that the corresponding solution stuck to some 
regularity island during its evolution. However this sticking is not permanent 
and the background diffusion grows when longer time spans are considered. 
The background values are shown in Fig. [6ja) only for the time span 10 
years to allow a better comparison of the numerical results with the dashed 
lines representing results from Chirikov's model, 

2. For Am D ~ 2, the background diffusion shows a discontinuity which, for the 
longest runs, reaches up to 14 orders of magnitude. This means that the cen- 
ter of the resonance domain becomes regular and the stochasticity remains 
confined to layers around the separatrix. This is the domain where Chirikov's 
slow diffusion theories are valid and where the results may be compared to 
the theoretical results obtained in Sect. 12.41 The integration time is a crucial 
factor in the detection of the slow diffusion. For instance, one may see that for 
simulations over only 10 years, the diffusion near separatrix is equal to the 
background diffusion for values of A mD close to 1, while for simulations over 
10 10 years, the equality is reached only for \m D = 9. 

3. In the last section of the Figs. [6] the solutions starting close to the separatrices 
show a diffusion equal to the background diffusion. The interpretations is that 
the stochastic layer is this case is so thin that the used initial conditions are 
no longer within them. (For that sake, the locus of the separatrices should 
be computed with very large precision. See e.g. Froeschle et al. 2006). One 
striking feature in this section is that an increase in the time span by a factor 
10 means a decrease in the background diffusion by a factor 10 3 . This is a clue 
for the fact that the solutions are dominated by periodic terms. Indeed, if we 
consider one periodic term with amplitude proportional to e and frequency u, 
its contribution to the average momentum in an interval [a, ft] is proportional 
to 



where At = ft — a. This integral is elementary and the integration of the result 
over all frequencies below a upper limit wi; m , gives 



b 




a 



— [si (6oj lim ) - si (ow Um )] 
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Fig. 6 Asymptotic behavior of the diffusion coefficient D3 for initial conditions over the 
separatrix and on the central part of the guiding resonance for (a) e = 0.05 and (b) e = 0.20. 
The dashed lines show the behavior predicted with Chirikov's theory. 



where si is the sine-integral function. The diffusion coefficients are given by 
the square of the average variation of the momentum divided by the total time 
(see Eqn. O and then D ~ At~ 3 . We also have D cc e 2 ~ J?" 2 ~ A~ 4 . The 
inclination -4 of the straight lines in the log-log plots can be easily checked. 

The diffusion of the solutions in the neighborhood of the separatrix may be deter- 
mined from Eqn. 1271 This equation involves the intensity of the perturbation (re- 
lated to A mD ) and two unknown parameters: the factor of reduction R mD and the 
factor of odd perturbations Rt- The factor of reduction corresponds to Chirikov's 
hypothesis of reduced stochasticity (due to holes, the solution does not fill the strip 
around the separatrix); the other factor comes from the fact that the perturbation 
is not even and thus the values of the diffusion coefficient are not the same for 
solutions in both separatrices (the solution may remain circulating near one of the 
separatrices at time different of the time it remain near the other). 

The results obtained with Chirikov are shown in Figs. [6] by dashed lines. In 
Figs.[6ja) tree different solutions are shown (calculated with the reduction factors 
indicated in the figure). The better agreement is obtained with R mr > = 0.25. The 
two values used for Rj- (0.6 and 0.9) give almost the same result, showing that 
the motion near deviation for the weakest perturbation (larger A mr > ) ■ In the other 
two figures, only the two solutions with R mr > — 0.25 are shown. 
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7 Conclusion 

Chirikov's theories provide heuristic tools to understand the diffusion observed in 
both eccentricity and semi-major axis of asteroids inside the (5, —2, —2) resonance. 
The multi-dimensional Hamiltonians of the three-body (three orbit) mean-motion 
resonances may be studied with the theories developed by Chirikov and collabo- 
rators, mainly because of the particular geometry of those resonances in the plane 
(a,e). The results obtained in this paper for the (5,-2,-2) three-body mean- 
motion resonance confirms the role of the resonances in the raising of diffusion 
across and along the main resonance as foreseen in Chirikov's theories. 

The diffusion calculations presented in this paper show that diffusion in semi- 
major axis is related with the diffusion in the across actions (pi,P2) while the 
diffusion in eccentricity is related with the diffusion in the along actions (p3,P4). 
The diffusion coefficient for the semi-major axis tends to small values showing 
that the variation of the semi-major axis remains small. It indicates the existence 
of barriers on both sides of the stochastic layer limiting the motion across the 
resonance. 

The comparison between simplified and complete model results shown that the 
diffusion in eccentricity is presumably due to the presence of at least one resonance 
driving the motion along the guiding resonance. This behavior is similar to the 
expected behavior of the Arnold diffusion, but, differently of it, the diffusion here 
is well apparent and the diffusion coefficients remain high. For this reason it was 
sometimes called Fast Arnold Diffusion (Chirikov and Vecheslavov 1989, 1993). 

The structure of the (5, —2, —2) resonance is formed by several overlapping res- 
onances, particularly at low eccentricities. Thus, diffusion across (5,-2,-2) reso- 
nance may be no longer limited to the thin chaotic layers (stochastic layers), but 
it fills the whole resonance zone. The diffusion along the resonance could be due 
to a multiplet. In this scenario, we have a random motion across the resonance, 
due to the overlap of several resonances belonging to a multiplet, and another, 
likely due to weaker resonances, which drive the diffusion along the guiding reso- 
nance. Arnold diffusion might occur inside the stochastic layer formed around the 
separatrix of the guiding resonance under the action of sufficiently weak pertur- 
bations. At variance, thick layer diffusion can appear for perturbation parameters 
in a broad interval. Although this mechanism show a similar exponential depen- 
dence of diffusion rate as a function of some system parameters, the mean rate of 
thick layer diffusion is generally larger than any theoretical estimation of Arnold 
diffusion. Therefore, it seems that for the real problem would be more appropriate 
to call the diffusion with another name - asymptotic diffusion of Chirikov- Arnold 
- due to the fact that this keeps some features of the Arnold diffusion, but late 
very well characterized by Chirikov like some distinct. 

As we mentioned before, we believe that the connection between rigorous in- 
vestigations concerning strictly Arnold diffusion and that observed in real physi- 
cal systems like this, is still an open subject. As Lochak (1999) pointed out, the 
global instability properties of near-integrable Hamiltonian systems are far from 
well-understood. It could almost be said that little progress has been made after 
pioneering work Arnold, and new ideas are definitely called for. 

Finally, the good results showed that the Chirikov slow diffusion theory can be 
used in broader investigations considering more resonances for (5, —2, —2), as well 
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applied for the others three-body (three orbit) mean motion resonances and also 
can include the inclination of the asteroid orbit. 



A Estimate of total variation of the momenta Pfc's 

To estimate the integral in (1201 we use the approach done by Chirikov (1979). In fact, the 
unperturbed separatrix is defined by 

ij>l x (t) = 4arctan (e ±f? G(t-t )^ _ ( 50 ) 

pf = ±2\M G \n G sm^, (51) 



where 

V G 



(52) 



' M G 

is the proper frequency of the pendulum Hamiltonian Hi . The double sign indicates the two 
separatrix branches: The positive sign correspond to the upper separatrix (0 < ipl x < 27r) , and 
the negative corresponds to the lower separatrix (— 2tt < ipf x < 0). The stable equilibrium 
points lie at ipi = ±7r, respectively. (This non usual separation of the intervals where the 
two branches are considered allows ipi and pi to have the same signal in each separatrix and 
simplifies the next calculations. For the usual presentation, the reader is referred to the study 
of the motions near the separatrix of pendulum in the Appendix B of Ferraz-Mello, 2007.) 
We introduce a time variable change r = Qq (t — t ), with tjif (t°) = i/>J = ±7r. Then, to 
> (lower separatrix), we have 

sin (t) = sin (£ m V>r (r) + A m r + ^) , (53) 

where i^^, = ^mV'i + w m (° + /9 m , with i/>i = i", and 

A m = ^- (54) 

Or, after expansion of the right-hand side, 

sin^ (t) = sin (£ m ipf x (r) + A m r)cos^ + cos(f m ^£ a! (r) + A m r) sin (55) 
When (1 55 It is substituted into H20II . the first term does not give contribution since, by symmetry, 

sin {UVi X (r) + A m r] dr = 0. (56) 



/ 



The contribution of the second term of d 5 5 II is determined by the relative signs of £m and A m . 
Using the absolute values to £ m and A m , we can introduce the Melnikov integral in the form 



/ 



+ °° 1 

cos(|£ m |V>rM± \\ m \r)dr = — A 2 | 5m | (T|A m |), (57) 



where j42|5 m | i s the Melnikov integral with argument ± |A m |. Then, the integral in J 20 1 ) i 



is 



r 



sinvm (*) dt = 77- sin Vm^2|e m | (=F |A m |) ■ (58) 



In the other branch of the separatrix, if)\ x has the signal changed, but the particular 
symmetry of this equation makes it invariant to the sign change of ipl x and, thus, one obtains 
the same result 1581 1. Indeed, if ipf x < the parity of cosine makes the integral H57H to be 



/ 



cos m x Ml T |A m | r) dr = — A 2 | 5m | (± |A m |) , (59) 
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where we used ipf x (r) = — \ipf x (t)|. Then, the variations in the actions pj. can be obtained 
introducing the result J59H into (1 20 P : 

4p fc (t)«-f- J2 Mm)V^sin^[A 2 | Cm |(|A ra |) + A 2Km |MA m |)]. (60) 

Chirikov (1979) estimated the diffusion using the result of the last equation. In order to simplify 
the theoretical estimate of the diffusion, Chirikov considered only even perturbing resonances 
and neglected the contribution of the perturbation with negative argument under the condition 
|A m | > 1. 

In the case of the three-body mean-motion resonance model, the perturbation are non 
even and it is not possible to neglect the contribution of perturbations for which |A m | is small. 
Then, each perturbation contributes differently when the motion lies close to a separatrix 
where A m > or A m < 0. Moreover, the odd perturbations in the Nesvorny-Morbidelli model 
makes necessary to take into account that the times of permanence of the motion near each 
separatrix are not equal. This is done by considering that the solution lies only a fraction of 
total time near the separatrix with A m > 0. To take into account this asymmetry we introduce 
the factor Rrp 

Rt = ^, (61) 

where T\ is the time that the solution stay in the neighborhood of separatrix with A m > 
and T is the total time. Then, the Eqn. H60II is rewritten as 

Ap k ~ 77- ^fe( m )QmSin^, (62) 

with 

Qm = V£ [R T A 2liml (|A m |) + (1 - R T ) A 2 | em | (- |A m |)] . (63) 

Equation (162 1 ) is valid for non-even perturbation and for small A m . In order to obtain estimation 
of H62B in terms of ordinary functions we must know the values of |A m | and |£ m |- In general, 
the relations to ^2|£ m | depend on the exponential term with argument |A m | (see Appendix A 
in Chirikov 1979). 
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